Relaxation times of kinetically constrained spin models with glassy dynamics 
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We analyze the density and size dependence of the relaxation time r for kinetically constrained 
spin systems. These have been proposed as models for strong or fragile glasses and for systems 
undergoing jamming transitions. For the one (FAlf) or two (FA2f) spin facilitated Fredrickson- 
Andersen model at any density p < 1 and for the Knight model below the critical density at which 
the glass transition occurs, we show that the persistence and the spin-spin time auto-correlation 
functions decay exponentially. This excludes the stretched exponential relaxation which was derived 
by numerical simulations. For FA2f in d > 2, we also prove a super-Arrhenius scaling of the form 
exp(l/(l-p)) < r < exp(l/(l-p)2). For FAlf in d=l,2 we rigorously prove the power law scalings 
recently derived in "|l3ll while in d > 3 we obtain upper and lower bounds consistent with findings 
therein. Our results are based on a novel multi-scale approach which allows to analyze r in presence 
of kinetic constraints and to connect time-scales and dynamical heterogeneities. The techniques are 
flexible enough to allow a variety of constraints and can also be applied to conservative stochastic 
lattice gases in presence of kinetic constraints. 



Despite a great deal of theoretical and experimental 
efforts, the main issues in understanding liquid-glass or 
jamming transitions remain unsolved. These phenomena 
occur in different systems: supercooled liquids, colloidal 
suspensions, vibrated granular materials 0, Char- 
acteristic features of a glassy behavior include dramatic 
slowing down of dynamics when a proper external param- 
eter is tuned (e.g. temperature for liquids) and the occur- 
rence of a complex, non exponential and spatially hetero- 
geneous Q relaxation process. Experiments show that, 
if an ideal glass transition occurs, it exhibits mixed first 
and second order features. Indeed, even if the divergence 
of the relaxation time signals a second order like tran- 
sition, no static diverging correlation length is detected. 
Furthermore relaxation times for fragile liquids diverge 
in a super-Arrhenius way: t ~ exp {A{T — Tc)/T) with 
A{x) I 00 as a; I 0, a signal of a cooperative relaxation 
on increasingly large scales as T decreases. 

A class of microscopic models that have been proposed 
in the attempt of understanding these phenomena are the 
so-called kinetically constrained models (KCMWijJaJa, 

S S la, 10, ii, 12, ii, M, i5, d6, 17, ,18, 1% m Mm 

and references therein). KCM are systems of particles 
on a discrete lattice with no static interactions beyond 
hard core and evolution is given by a Markovian stochas- 
tic dynamics. The elementary move is either a jump 
of a particle to a nearby empty site in the conservative 
(Kawasaki) case or the creation/destruction of a parti- 
cle in the non conservative (Glauber) case. In both sit- 
uations the associated rates, if positive, verify detailed 
balance with respect to a simple product measure. How- 
ever, and this is the most interesting feature, a move 
at site X can occur only if the configuration around x 
satisfies certain constraints. The latter mimic the local 
constraints that may occur in the physical systems and 



which may cause the dynamical arrest. Indeed, for a 
proper choice of the constraints, KCM show glassy fea- 
tures including stretched exponential relaxation, super- 
Arrhenius slowing down and dynamical heterogeneities. 
Several works and a great deal of numerical simulations 
have recently been devoted to understand the mechanism 
inducing these glassy properties and to derive the typical 
time-scales as well as the asymptotic form of correlation 

functions (iJl^ H 11 E [H III 111 III III III III 

ilSi. il^. I20I I211 1231'). Numerical simulations are however 
very delicate because of the rapid divergence of r as the 
density is increased. Furthermore, finite size effects often 
have non trivial scalings. In this work we present a new 
general probabilistic technique through which we obtain 
rigorous results on the dependence of t on the size of the 
lattice and on the particle density, p. As a by-product we 
also obtain meaningful bounds on the long-time behavior 
of time correlations and persistence function. 

We first introduce the models and our results and then 
sketch the technique, referring to 23] for rigorous proofs. 
We focus on non conservative models, also known as ki- 
netically constrained spin models (KCSM) or facilitated 
spin models. Conservative models will be analyzed else- 
where |2^. KCSM are defined as follows. Each site x 
can be occupied, rj^—l, or empty, rjx^^O, and it changes 
its current state with rate [(1 — p)rix -t- p(l — r]^)] fxiv)- 
Here /a, (77) does not depend on rj^ and it encodes the ki- 
netic constraint. Thus detailed balance is satisfied w.r.t. 
the Bernoulli product measure, ^p, at density p. KCSM 
can be divided into two classes: (i) non-cooperative and 
(ii) cooperative. In a non-cooperative model it is possible 
to completely empty any configuration containing some- 
where a finite seed of vacancies. In a non-cooperative one 
that is not the case. The class (ii) can be further divided 
into: (iia) models that are ergodic in the thermodynamic 



limit at any p < 1; (iib) models that display an ergodicity 
breaking transition at a critical density pc < 1- For the 
latter, above pc, there exists an infinite spanning cluster 
of particles mutually blocked by the constraints. Finally 
(iib) can be classified as discontinuous or continuous ac- 
cording to the character of the percolation transition of 
the blocked structure. 

Among non cooperative models we consider the one 
spin facilitated Fredrickson- Andersen model (FAlf), 
which recently received a renewed attention as a model 
for strong glasses ^ move is allowed only if at 

least one of the nearest neighbors is empty: f'n(Ti)=l if 
T,y n.n.xi^ " Vx) > 0, fx{v)=0 Otherwise. In [uj a dy- 
namical field theory was derived which gives {q=l—p): 
T (xl/q'' with z=3 for d^l, z^2 + e{d) with e(2) ~ 0.3, 
e(3) ~ 0.1 and e{d > 4) = 0. A recent exact mapping into 
a diffusion limited aggregation model, its renormalization 
and a careful treatment of the symmetries involved 
gives instead d—2 as the upper critical dimension and 
e(d)=0 in d > 2. 

Our results for FAlf are the foUowings. If t denotes 
the inverse of the spectral gap of the Liouvillian operator 
generating the stochastic dynamics, i.e. the worst relax- 
ation time on all one time quantities, we get: r cx l/q^ 
in d = 1, < T < l/(g^logl/g) in d = 2 and 

l/g-'^+^Z'' < T < 1/q^ in d = 3. These rigorous results 
lead to e(2) = and e(3) < 0, disproving the findings in 
[m and confirming those in 0|. Moreover our method 
allows us to identify explicitly the slowest modes which 
dominate relaxation. Furthermore, it is easily adapted to 
any possible choice of non cooperative constraints, e.g. to 
models in which more than one spin is needed to empty 
the whole lattice. We can also treat models with asym- 
metric (e.g. East model ) or partially asymmetric con- 
straints (J7|) which have been proposed to model fragile 
glasses PJJ and the strong/fragile crossover. In partic- 
ular for East, we obtain logr ~ {\og{l/q))'^ as in 
Finally we can deal with the persistence function 

Fit) = J dpp{r,m P[77o(5) = 77o(0) < t] . 

In great generality |2^ we prove that whenever the global 
relaxation time r is finite, F(t) decays exponentially fast. 
In particular, for FAlf we obtain Tp < q~^T for small q, 
where Tp is defined via F{tf) = e^^. Although these 
findings disprove any stretched exponential decay at large 
times, for FAlf in d = 1 our bound Tp < 1/q* does not 
preclude the stretched form [ill F{t) ~ exp(— y^^^) 
when q I and q'^t ~ 0(1). This is indeed the regime 
where numerical simulations are performed [l^ . 

Among cooperative models without transition (iia) we 
consider the two spin facilitated FA model (FA2f) on 
square and cubic lattices. Here the constraint requires 
that at least two of the surrounding sites are empty in 
order for the rate to be non zero. Originally ^ a dynami- 
cal transition at a finite density was predicted but later it 



was shown that it cannot occur Q,^,!!^. Still, FA2f is 
relevant for fragile glasses since both super- Arrhenius di- 
vergence for r and stretched exponential relaxation have 
been detected |l[ll[lllil[llElIIi- Indeed, both the 
spin-spin time autocorrelation and the persistence func- 
tion are fitted with exp (— (t/r)^) with P decreasing as 
the density p is increased 0, 0, . As pointed out 
in Q , beyond the general recognition that FA2f behaves 
like fragile glasses, little is known about the scaling of 
r. Different fits have been proposed: Adam-Gibbs form 
[13, Vogel-Fulcher form [3 and exp(c/g) [3. The lat- 
ter, which corresponds to a super- Arrhenius form when q 
is rewritten in terms of the temperature, is supported by 
the conjecture that relaxation occurs via the diffusion of 
critical droplets of size 1/q over distances exp(c/g) [l^ . 

Our results for FA2f in d = 2 and d = 3 are the foUow- 
ings. We prove that exp{l/q) < r < exp(l/g^) which 
establishes a super-Arrhenius scaling and excludes the 
Vogel Fulcher form Furthermore, since our results 
hold uniformly on the system size, we obtain that relax- 
ation in infinite volume is purely exponential at any fixed 
q. Furthermore, as for FAlf, we get a strict bound on the 
crossover regime where a stretched form may occur if we 
let simultaneously q i and t 1 oo. These findings con- 
tradict the asymptotic stretched exponential relaxation 
in 0, 0, 0| . The fit with a stretching exponent /3 < 1 
should be due to the rapid divergence of t, which was 
also a fitting parameter. 

Among (iib) models with continuous transition, we 
consider the two dimensional North East model (NE). 
For NE both the up and right neighbors should be empty 
in order for a move to be allowed and a continuous transi- 
tion occurs at the critical density of directed percolation, 
Pdp [2(^1 . Finally, for (iib) with a discontinuous transition, 
we consider the recently proposed Knight model (2^ . We 
refer to [l^ for the precise definition of the constraints; 
we only recall that the transition occurs again at p^p and 
has the remarkable features of an idea7 glass transition. 
Indeed a finite fraction of the system is frozen at pdp and 
the size ^ of blocked structures diverges faster than any 
power law as p t Pdp- This leads to a super-Arrhenius 
divergence if r ~ is assumed. 

Our main results for NE and Knights are the foUow- 
ings. For p < Pdp we identify a constant which bounds t 
uniformly on the system size, L. Therefore, even if an er- 
godicity breaking transition occurs, relaxation is purely 
exponential in the ergodic region. Moreover, we get that 
T diverges when p y pdp and r oc exp(L c(p)) for p > pdp- 
This result provides a possible test on the validity of 
Knights for systems undergoing jamming transitions. 

Let us now sketch the main ideas beyond our approach, 
rigorous proofs will be reported elsewhere [i^. Here we 
focus on d = 2. Results in d = 1 can be derived analo- 
gously; we will comment at the end on d = 3. We first 
introduce a new model referred to as the General Model 
(GM) and derive results for its relaxation time, tgm- 
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FIG. 1: a) Block dynamics for GM: percolating path of good 
sites (•) required to renew configuration on 3. b) Blocking 
event for FA2f. • (o) stand for particles which do (do not) 
belong to the backbone. Sites inside dotted line form one of 
the sequences of > 5Lc/2 sites to be emptied before O. 



Then we show how to map a given KCSM into GM by 
a renormahzation technique. GM can be described as 
follows. Instead of two- valued occupation variables con- 
sider N- valued variables, rix € S = (0, . . .TV — 1), and 
a probability measure on S. We identify a subset G 
of S which we call the good event and we declare good 
a site X if Ux ^ G. The dynamics is defined as follows. 
Each site x waits a mean one exponential time and then 
Ux is refreshed by a new value n'^ sampled from ly, pro- 
vided its three North, North-East and East neighbors 
(i.e. X + ei, X + ei + 62, X + 62, being the unit basis 
vectors) are good. If this constraint is not satisfied, 
remains unchanged. We consider GM on a square lat- 
tice of linear size L with good boundary conditions 
on the top and right boundaries to ensure ergodicity. In 
order to evaluate tgm (L) we follow the bisection method 
[2^. Partition Al into four blocks as in fig^) and de- 
fine the following auxiliary accelerated dynamics. Each 
block waits a mean one exponential time and then its 
configuration is replaced by a new one chosen according 
to the product equilibrium probability given by v. On 
the top right block (block 2 in fig^)) this move is al- 
ways allowed because of the boundary conditions. For the 
others, a constraint should be satisfied: on an 1-shaped 
frame of width , d < 1, there should be a percolat- 
ing cluster of good sites in the current configuration as 
in fig^Jt). In other words, the constraint requires the 
good GM boundary conditions on the blocks 1,3,4 (on 2 
they are automatically guaranteed). The relaxation time 
tqm{L) is then bounded by the product 

tgm[L) < rbiock(^) TGM (i/2) 

with Th]nr\<(L) the relaxation time for the block dynam- 
ics |2^. The above inequality corresponds intuitively to 
a two steps relaxation: first on the block scale, then in- 
side each individual block. Notice that, trivially, tgm{L) 
is smaller than the worst case over the boundary condi- 
tions of the relaxation time on scale L/2. However, such 
inequality is useless because without the good b.c. GM 
dynamics is not ergodic (i.e. r is infinite). In other words, 
working with a constrained block dynamics is essential. 
Then, by dividing Ai/2 into four blocks and so on up to 



where c is a finite constant and the product contains 
O(logL) terms. Therefore we get a bound for tgm{L) 
which does not dependent on L provided that the prod- 
uct over the Tbiock's converges. In turn this occurs if the 
probability that a site is good, p, is larger than a finite 
threshold, pc < 1- Indeed one can easily show that, for 
p ~ 1, rbiock(i) — (1 -|- exp(— mL'^)) where exp(— mL'') 
comes from the probability that the constraint for the 
block dynamics on scale L is violated. 

Let us now explain how to apply such a result to the 
KCSM models described above. We consider for definite- 
ness FA2f in d = 2 with empty boundary conditions on 
the top and right borders. We say that a region V is inter- 
nally spanned for the configuration rj if, when consider- 
ing all the sites outside V occupied, rj can be completely 
emptied using (internal) allowed moves. The probabil- 
ity for Al to be internally spanned has been evaluated 
in the context of bootstrap percolation: it goes to one 
exponentially fast when L exceeds the crossover length 
Lc — exp{c/q) with c — 7r^/18 26]. Divide now the 
square lattice into disjoint blocks of size kLc, A: 3> 1, 
and identify each block by the coordinates of its cen- 
ter. On each of these renormalized sites we define a 
configuration space S = {0,1}^', with Nb the number 
of sites in a block. In S we identify the good event G 
as the set of configurations such that the block is inter- 
nally spanned. Because of our choice of k the probability 
that a renormalized site is good is ~ 1. We can thus 
run the GM dynamics on the renormalized lattice and 
get sup^, tqm{L) ~ 1. Since in GM dynamics blocks 
are updated only if "enough" surrounding are internally 
spanned, by the two steps relaxation argument we get 

t{L) < TGM{L)T{kL,) - r(fcL,). 

where t{L) is the relaxation time for FA2f. In the infinite 
volume limit and for any p < I this leads to an exponen- 
tial relaxation for all one time functions as well as for 
the persistence function using a Feyman-Kac bound . 
At the same time the density dependence of t is com- 
pletely encoded in the relaxation time on a lattice of size 
Lc{p)- To evaluate the latter we reduce the scale from 
Lc to l/g^ via a strategy similar to the previous one. 
However, on scales smaller than Lc the event "the block 
is internally spanned" becomes very unlikely and we are 
forced to make a different choice for the constraints of the 
renormalized dynamics in order to keep tgm — 1- Our 
choice for the good event G is suggested by the follow- 
ing two observations: (i) any empty segment of length £ 
can be displaced by one step in a given direction if there 
is at least one vacancy on the adjacent segment in that 
direction; (ii) the probability that there exists a fully oc- 
cupied segment of length £ inside a critical square of size 



Lc is very small as soon as I ^ l/?^- Thus, we choose 
good events which force the not fully occupied condition 
on segments of length l/q^ . By applying a bisection pro- 
cedure analogous to the one used for GM together with 
elementary paths arguments, we get 

r(L,) < cL,T{l/q^) 

Finally we bound T{l/q^) with the highest entropy cost 
to get tfa[L) < cLcexp(l/g2) = 0(exp(l/g2)). 

The results for the other 2-dimensional models (FAlf, 
NE and Knight) are obtained analogously. Furthermore, 
for a generic choice of the constraints, we find that t is 
dominated by t{Lc) where Lc is the crossover length for 
a region to be internally spanned (with the chosen con- 
straints). Therefore Lc, which can be determined by a de- 
terministic bootstrap-like procedure, is the relevant size 
over which numerical simulations should be performed. 
By analyzing the constraints one can also further reduce 
the problem to a much smaller scale, as we do for FA2f. 

Lower bounds for r, directly in the thermodynamic 
limit, can be established either via a suitable choice of 
test functions in the variational characterization of r 
or using our general result for the persistence function: 
F{t) < exp{—tq/T) ■ Consider an event B, called the 
blocking event, and let Pb {t) be the probability that the 
origin is occupied for all times s < t, worst case among 
all the starting configuration in B. The inequality for 
F{t) implies fXp{B)PB{t) < exp(— ig/r). We define the 
blocking event B as the set of configurations for which, 
after standard bootstrap percolation (i.e. all possible 
FA2f moves in sequence) inside A^^^, a backbone of par- 
ticles containing the origin survives. By choosing 6^1 
we have ^ip{B) ~ 1. In infinite volume this backbone will 
eventually get unblocked thanks to the vacancies outside 
Aji^. However, this requires an ordered sequence of at 
least 5Lc/2 moves (figIT]D). Thus, Psit = eSLc) ~ 1 for 
sufficiently small e. Therefore 0(1) < exp(— ig/r) for 
t ~ eSLc, i.e. t > 0{Lc)- In rf = 3 lower bounds can 
also be obtained as before: using the results in for 
Lc we get T > exp{c/q) for FA2f and r > expexp(c/g) 
for FA3f. The latter is in agreement with the conjecture 
in 1 141 and disproves the possibility of an Adam Gibbs 
form [13 which corresponds to r ~ exp{c/q). Finally 
the sharp lower bounds mentioned before for FAlf are 
obtained by a careful choice of a test function. 
In summary, we have developed a technique which 
allows to obtain rigorous results on r for KCSM via 
the knowledge of the typical regions which have to be 
rearranged to perform a movement. In other words, we 
have been able to connect rigorously time-scales and 
dynamical heterogeneities, a subject which have recently 
received great attention in the glass community [2|. 
The main new results established via this technique are 
super- Arrhenius behavior and exponential relaxation 



for cooperative models. Moreover our method suggests 
the typical length scale on which numerical simulations 
should be carried on. As for future developments, we 
believe that our technique can be applied to finite 
dimensional models other than KCM, e.g. models with 
static interactions beyond hard core and to conservative 
KCM models. 
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